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ABSTRACT 


The goals of this thesis are 1) to design and test a two-axis surface wave source 
and two-axis surface wave receiver and 2) investigate surface waves to detect buried 
objects in water saturated sand. Results of measurements confirm the ability to generate 
particle motions in water saturated sand consistent with surface wave excitation. 
However, limitations in the size of the test tank prohibit a thorough investigation of the 
ability of the source and receiver to selectively excite and detect surface waves 
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I. INTRODUCTION 


A. OBJECTIVE 

The purpose of this research was to design, build, and test a two-axis source and 
two-axis receiver to investigate the use of surface waves to detect buried objects in water 
saturated sand. 

B. MOTIVATION 

Previous research at the Naval Postgraduate School evaluated the feasibility of a 
phased array surface wave source to detect buried land/sea mines [Ref 1]. The source 
and geophone combination used in that work were single-axis, constrained to vertical 
motions. That work demonstrated target localization using surface wave scattering, and 
beamforming with phased array techniques; however, the transmitted and received waves 
were not folly characterized. A logical progression of that research was to use a two-axis 
source and two-axis geophone combination to specifically excite and acquire surface 
waves. This two-axis source and two-axis geophone combination could then be used to 
investigate the feasibility of using elementary multi-axis signal analysis to enhance our 
understanding of surface wave propagation in water saturated sand. The purpose of this 
thesis was to build upon the results of the aforementioned research with a two-axis source 
and two-axis geophone combination. Similar techniques have been used by others [Ref 
2 ]. 
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C. THESIS OUTLINE 

The following topics will be discussed in this thesis: 

1. Description of the test environment and data acquisition system. 

2. Design and calibration of the two-axis geophone. 

3. Design and characterization of the two-axis source. 

4. Comparison the single-axis source/two-axis geophone combination to the two- 
axis source/two-axis geophone combination. 

5. Investigation of the feasibility of using polarization filtering and Cepstrum 
analysis with the single-axis source/two-axis geophone combination and the 
two-axis source/two-axis geophone combinations to detect buried objects in 
water saturated sand. 
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11. APPARATUS 


The purpose of this chapter is to 1) describe the test tank, 1) describe the 
LAB VIEW data acquisition system, 3) describe the design and calibration the two-axis 
geophone, and 4) describe the design of the two-axis source. 


A. TEST TANK 

All of the experiments associated with this thesis were conducted in a round, 1.8 
radius, 1.2 m deep redwood tank. The tank’s redwood bottom covers a hard cement 
laboratory floor. This tank was filled with approximately 17 tons of medium grain (0.6 
mm) beach sand from Moss Landing in Monterey County, California [Ref 1]. There was 
a 3.7mx5cmx3.2cm aluminum beam mounted 10 cm above the sand. This beam was 
used to mount the source above the sand (see Fig. 2.1). The tank was conditioned prior 
to each experiment. This conditioning included uniformly wetting the surface with a 
garden hose with spray attachment, and uniformly pounding the surface of the sand with a 
50 cm X 50 cm square board attached to a wooden pole to ensure uniform packing of the 
grains near the surface. 
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Figure 2.1 Diagram showing top view of the water saturated sand-filled redwood test 
tank. The two lines leading to/ffom the source and receiver emphasize that both are 
two-axis devices. 
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B. DATA ACQUISITION SYSTEM 


Most of the data acquisition associated with this thesis was performed with 
LAB VIEW version 3.1.1. A product of National Instruments Corporation, it is a 
graphical, general-purpose programming language with extensive libraries of functions, an 
integral compiler and debugger, and an application builder for stand-alone applications 
[Ref 3]. In addition to the software, three plug-in boards were utilized; the NB-MIO- 
16XL Data Acquisition board, the NB-AO-6 Analog Output board and the NB-DMA- 
2800 Direct Memory Access board. The LABVIEW software was used on a Power 
Computing model 100 computer (see Appendix A). Source code for post-acquisition 
processing was written using MATLAB 4.2. 

C. TWO-AXIS GEOPHONE DESIGN PARAMETERS 

As in Ref 1, the design requirements for the two-axis geophone to detect 
vibrations in the sand were the following: 

1. It should have adequate sensitivity to ground motion generated by the source 
within the ranges provided by the tank. 

2. It should have adequate sensitivity in the frequency band of interest, nominally 
200 to 400 Hz. 

3. It should have good coupUng to the sand that could give repeatable results. 
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Additional requirements for this thesis were; 


4. It should be capable of sensing both vertical and horizontal components of the 
signals. 

Figure 2.2 shows a diagram of the two-axis geophone. The sensors are Applied 
Magnetics Geo Space model number HS-1 miniature geophones (see data sheet in 
Appendk B). The two-axis geophone housing is made of machined aluminum. The 
vertical geophone is secured to the housing by a threaded rod manufactured into the 
geophone case. The horizontal geophone is secured by set screws in the aluminum 
housing. The 5.2 cm sand ‘coupling peg’ was the same coupler used on the model 
S4408A sensors used in Ref 1. 


D. TWO-AXIS GEOPHONE CALIBRATION 

The two axis geophone was calibrated using a HP3562A Dynamic Signal 
Analyzer, (APS) Acoustic Power Systems Inc. Dual Power Amplifier, APS Model 1201 
Perma-D 5 me Shaker Table, and a calibrated model J353B03 PCB Piezotronics 
Accelerometer. The horizontal axis geophone was calibrated by mounting the shaker 
horizontally on a stable base. The geophone and accelerometer were securely mounted on 
the shaker table with a specially manufactured mounting plate. They were mounted with 
their sensitive axis aligned horizontal. The swept sine output fi'om the HP 3562A was 
amplified by the APS Dual Power amplifier and applied to the shaker table. The frequency 
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Figure 2.2. Diagram showing the design of the two-axis geophone. 
The hatched portions indicate the individual single-axis geophones 
mounted in an (unhatched) aluminum block. 


7 








range of the sweep was 0-1000 Hz. The resultant accelerometer and geophone responses 
were applied to channels one and two of the dynamic signal analyzer. The built-in 
frequency response was used to determine the frequency response from these signals. Data 
was read from the display, using cursors, in 20 Hz increments and entered into a Matlab 
file. The calibration curve shown in Fig. 2.3 was computed in the following manner. The 
velocity of the accelerometer must equal that of the geophone 

Ua=Ug. (2.1) 

The velocity is related to acceleration by u = a/©. The sensitivity of the accelerometer Sa 
is given in mV/g, where g = 9.81 m/s^. The voltage output of the accelerometer Va is 
converted to a velocity according to 

Ua= (9.8 m/s^/g) (Va/Sa©), (2.2) 

where V* is in mV. The voltage output of the geophone is converted to velocity 
according to 


% Vg/Sg 


(2.3) 
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where Vg is in mV and Sg is the sensitivity of the geophone in mV/m/s. Setting the two 
velocities equal and solving for Sg gives 

Sg=(VgA^a)(a)/9.81m/s^/g)Sa. (2.4) 

The vertical axis geophone was calibrated in the same manner except the shaker was 
mounted vertically. 

As can be seen from Figs. 2.3a and b the vertical and horizontal calibration curves 
are similar over our anticipated operating frequency range. Therefore for comparison 
purposes, equal voltages from the two geophones correspond to approximately equal 
velocities. In the results presented later, the output voltages will be used rather than 
velocities. 

E. TWO-AXIS SOURCE DESIGN PARAMETERS 

The two-axis source design was guided by some of the same requirements as in 

Ref 1: 

1. It should generate signals with sufficient amplitude to allow detection of the 
surface waves anywhere in the test tank. 

2. It should provide a repeatable waveform. 

3. It should be capable of a wide range of repetition rates. 

4. It should have good ground coupling that would last for an extended period. 
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5. It should provide clean signals over a range of frequencies from 
approximately 200-400 Hz. 

The focus of this thesis required the following design criteria: 

6. It should be capable of generating two-axis motion. 

7. It should be portable and capable of being quickly moved and reset in order 

to facilitate conducting experiments of various source/geophone configurations 
in relatively short periods. 

Figure 2.4 shows the final design. Two four inch diameter model 4A9, high cone 
excursion, electrodynamic speakers were mounted at right angles to one another on an 
aluminum support plate (see data sheet in Appendix C). Connecting rods coupled the 
drivers to a 7.3 sq. cm source foot. The 7.5 cm rods were attached to the foot with 
machined holes and set screws. The mounting block’s dimensions were 2.5 cm x 1.0 cm x 
0.75 cm. The block was epoxied to the foot. The foot consisted of an array of 25,4 cm 
long nails mounted on a 2.7 cm x 2.7 cm brass plate and arranged in a square pattern, 
spaced by 0.6 cm to improve contact with the sand. The 7.5 cm rods were attached to the 
speaker cones by threaded aluminum caps that were epoxied to the speaker cones. 


11 



Figure 2.4. Diagram showing the design of the two-axis source. The two electrodynamic 
loud-speakers are mounted at a right angle. Two connector rods transfer the speaker’s 
diaphragm motion to the “bed of nails” foot. 
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F. SUPPORTING EQUIPMENT 


The typical experimental arrangement is shown schematically in Fig. 2.5. The two- 
axis source was driven by the outputs of two HP3314A function generators, each 
amplified by one channel of a Techron model 5507 power amplifier. One HP 3314A was 
triggered by the other, so that the relative phases of the outputs could be adjusted. 

The two geophone outputs were sent to bandpass filter-preamplifiers (one 
Stanford SR 560 and one Ithaco 1202). The outputs were sent to an oscilloscope and the 
A/D boards. 
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Figure 2.5. Block diagram showing the typical electrical equipment used in the 
measurements. The relative phases of the outputs of the two function generators can be 
adjusted over the range-180 to 180® 
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III. TWO-AXIS SOURCE CHARACTERIZATION 


The purpose of this chapter is to 1) describe the characterization of the two-axis 
source, 2) investigate frequency dependent attenuation in the test environment, and 3) use 
polarization filtering to fijrther investigate the ability of the two-axis source to excite 
surface waves. 

A. TWO-AXIS SOURCE CHARACTERIZATION 

Two PCB model 303AOS miniature accelerometers were mounted vertically and 
horizontally on the foot of the two-axis source. The source was driven with a 50 cycle, 
400 Hz pulse Avith pulse repetition frequency of 1 Hz. The phase relationship between the 
drive signals was varied and the resultant source-foot vibrations were sensed by the 
accelerometers and acquired with LAB VIEW. This experiment was repeated for two 
conditions; for the foot in air and for the foot coupled to the sand. Signals were also 
acquired from the two-axis geophone at a range of 20 cm from the source. The vertical 
signals were plotted versus the horizontal signals in order to compare the in-air 
accelerometer signals, the in-sand accelerometer signals and the received two-axis 
geophone signals. 

Figures 3.1 and 3.2 show a series of plots of the in-air and the in-sand 
accelerometer signals for different phase relations. The notation ’Rm5L90’ indicates that 
the right channel was driven by a signal with phase minus five degrees and the left channel 
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was driven by a signal with a phase of positive 90 degrees. The other plot labels can be 
interpreted in the same manner. These x-y patterns show that the two-axis source motion 
can be manipulated to produce a range of foot motions for various drive phases. The sand 
motion detected by the two-axis geophone is depicted in Fig. 3.3. The geophone signals 
were dominated by vertical motion. Although it was not possible to generate purely 
horizontal motion, it was possible to introduce a phase shift between vertical and 
horizontal motions. The reason for the relative lack of horizontal motion control is not 
clear. It may be due to restrictions of the foot motion by the sand or the connecting rods, 
or it may be due to some more fundamental limit. The foot is more free to move 
vertically. However, Fig. 3.3 shows that it is possible to generate near-elliptic motion in 
the sand. 

The motion of the sand as detected by the two-axis geophone is also dependent on 
the relative azmuthal location of the geophone with respect to the source. Figure 3.4 
shows the resultant patterns from signals received at 0, 90,180, and 270 degrees relative 
to the source at a range of 20 cm from the source. The 0° reference was a line 
perpendicular to the plane containing the connecting rods (see Fig. 3.5). As the geophone 
was moved from location to location the horizontal axis was pointed in a direction parallel 
to the connecting rod plane. AH subsequent experiments were conducted with the 
geophone in the 90® orientation. 
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Figure 3.4. Plots of the sand motion detected by the geophone for four different a 2 muthal 
orientations of the geophone relative to the source. The orientations are define in Fig. 3.5. 
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Figure 3.5. Top view showing how the azmuthal orientations used in Fig. 3.4 are defined. 
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B. RANGE DEPENDENCE OF SIGNALS 


Figure 3.6 shows two representative vertical and horizontal geophone signals 
measured at a range of 20 cm from the source. The two-axis source was driven with 
single-cycle, 90 degree phase shifted, 400 Hz signals with a 1 Hz repetition frequency. The 
data shown in Fig. 3.6 show a number of the features. First, the response in the sand was 
not a single cycle. Also the vertical and horizontal signals are very different. There is 
evidence of multiple arrivals, especially in the horizontal signal. The structure of the 
received signal is very range dependent. In an effort to demonstrate this, the waveform 
amplitude corresponding to the feature pointed out in Fig. 3.6, was plotted vs. range for 
signals in the frequency range 200 to 800 Hz and ranges 20 to 120 cm from the source. 
Figure 3.7 shows the results for frequencies 200, 400, and 600 Hz. The expected 
consistent monotonic decrease in amplitude with range was not observed. Instead, the 
irregular amplitude variations with range indicate that this tank does not model free field 
conditions. We decided to use 400 Hz for further investigation. 





0.02 0.03 0.04 0.05 0.06 

Time (seconds) 




Figure 3.6. Graphs of the output of the vertical (top) and horizontal (bottom) geophones 
as fiinctions of time. The source was driven in phase with a single 400 Hz cvcle 
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Figure 3.7. Plots of the amplitude of the vertical (top) and horizontal (bottom) geophones 
as functions of range for 200,400, and 600 Hz drive frequencies. 
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C. POLARIZATION FILTERING 

One way to detect surface waves in the presence of other modes of propagation 
(such as compressional and shear) is to use polarization filtering. As described in Ref 2 , 
the function of the polarization filter is to make use of the horizontal/vertical phase offsets 
characteristic of surface waves to preferentially extract surface wave signatures but reject 
signals due to other propagation modes. 

The first step in the analysis was to convert the real horizontal signal x(t) into a 
complex analytic signal X(t). X(t) is defined as 

X(t) = x(t)+jx'(t), (3.1) 

where x'(t) is the Hilbert transform of x(t) defined through 

x'(t) = l/7tJx(t)/(t-T)dT. (3.2) 

The vertical signal was similarly converted to an analytic signal Y(t). For monofrequency 
signals, the two functions X(t) and Y(t) should differ by a constant phase angle, 
characteristic of surface waves. The product X*(t)Y(t), were * denotes the complex 
conjugate, yields a complex term whose imaginary part is sensitive to the phase difference 
between the two signals. Im[X*(t)Y(t)] is zero if x(t) and y(t) are in phase; it is nonzero 
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otherwise. If the phase is such that the particle motion is prograde (retrograde) the 
imaginary part is positive (negative). 

To test this technique, this analysis was performed on two sets of computer 
generated waveforms. Figures 3.8 and 3.9 show the results. In the upper trace of Fig. 3.8, 
there are two waveforms (solid and dashed) that consists of two gaussian shaped pulses. 
The frequency of the signals in each of the pulses is 400 Hz. The relative phase between 
the first pulses of the two signals is 90®. The second pulses are in phase Avith each other. 
The lower trace shows the results when this analysis was performed on the signals in the 
upper trace. Note that the in-phase portions (the second pulses) of the signals were 
eliminated. The out of phase portions result in an envelope with period proportional to 
the period of the two out-of-phase (first pulses) signal portions. Figure 3.9 shows the 
results of polarization analysis on signals similar to the signals in Fig. 3.8. The only 
difference is the relative phase of the first pulses is 270° (-90°). Note that the 90° shift 
results in a positive output, while the -90° shift results in a negative output. 

When this algorithm was applied to a set of received horizontal and vertical signals 
as in Fig. 3.6 the results were very different from the theoretical results. This lead to a re¬ 
examination of the received signals. Upon further review, a difference was observed in the 
received vertical and horizontal signal’s frequency. This frequency difference is evident 
upon careful review of Fig. 3.6. Apparently the two-axis source has different frequency 
and phase responses in the vertical and horizontal directions. Because there is no 
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consistent phase relationship between frequencies of two different frequencies, 
polarization filtering cannot be used to detect surface waves. 

To demonstrate this problem, polarization analysis was performed on the two dual¬ 
envelope signals shown in Fig. 3.10. The upper trace shows two computer generated 
signals. The dashed signal consists of two gaussian envelopes. The frequency of the 
signal in the first envelope is 300 Hz and the frequency of the signal in the second 
envelope is 400 Hz. The solid line signal also consists of two gaussian shaped pulses. The 
frequency of the signals in each of the pulses was 400 Hz. The lower trace of Fig. 3.10 
shows the results of appl 5 dng this algorithm on these computer generated signals. One can 
see that polarization filtering gives a bipolar output for these signals. Figure 3.11 shows a 
comparison of the results of polarization analysis on the computer generated signals in Fig. 
3.10 and the real acquired signals. These signals are very similar. These results seemed to 
indicate that the polarization analysis conducted on the real signals revealed different 
frequencies in the received vertical and horizontal signals. 

The problem of two different frequency responses can be overcome by driving the 
source long enough to reach steady state. This is demonstrated in the next sequence of 
figures. Figure 3.12 shows the outputs of the vertical and horizontal accelerometers 
mounted on the foot for two different drive conditions, 1 cycle and 10 cycles. As might be 
expected, the accelerometer outputs for the 1 cycle drive roughly resemble decaying 
oscillations. However, the frequencies of the oscillations and the decay time were very 
different for the two axes. There is little similarity in the 1 cycle signals, whereas the 10 
cycle signals shows reasonably good similarity. Figure 3.13 shows the outputs of 
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Figure 3.10. Same as Fig. 3.8 except the signals in the first pulse (top) have ditterent 
fi-equencies rather than different phases. The output of the polarization filter (bottom) is 
bipolar in this case. 
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geophones for the same drive conditions as in Fig. 3.12. One notable difference is in the 1 
cycle data. The horizontal geophone signal lasts much longer than the foot motion. The 
reason for this difference was not clear. It might indicate either a reverberant 
environment, or resonance(s) in the source mounting beam, the source connecting rods, 
or the speakers. Figures 3.14 and 3.15 show the power spectral densities of the horizontal 
1 cycle and 10 cycle geophone signals. Figure 3.14 shows the vertical data. Fig. 3.15 the 
horizontal data. The main conclusion is that given a sufficiently long input signal, the 
source will produce the same frequencies in both axes. 

For signals of approximately 10 cycles or more, or approximately steady-state 
conditions, the response was closer to the excitation frequency (which was 400 Hz in this 
case). In short, to have a predictable transmission to the sand, the excitation signals’ pulse 
period must approximate steady state conditions in the source. However, assuming a 
propagation speed of 100 m/s, the wavelength at 400 Hz is 25 cm. Therefore, the length 
of a 10 cycle pulse is 2.5 m, which is comparable to the dimensions of the test tank. The 
conclusion to be drawn is that the test tank was too small to perform an adequate 
investigation of surface waves. 

These results were further complicated by the following; assuming an average 
group speed range between 50 m/sec and 100 m/sec and a minimum of 10 cycles per 
pulse, choosing the worst cases, the wavelengths of the source signals could range from 
2.5 to. 5 meters in space. These wavelengths are too large to resolve any features in the 
tank. Specifically, polarization analysis, or any algorithm that required frequency stability. 
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Figure 3.13. Graphs showing the output ot the geoptione tor me same conaiuons as m 
Fig. 3.12. 
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Figure 3.15. Same as Fig. 3.14 except for the horizontal geophone. 


Normalized Power Spectrum for Horizontal Geophone Signal (1 cycle) 



Normalized Power Spectrum for Horizontal Geophone Signal (10 cycles) 





















































could not be utilized with this two-axis design in this test tank. The same analysis was 
performed on signals excited by the single-axis source with similar results. 
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IV. SCATTERING EXPERIMENTS 

A. COMPARISON OF SINGLE-AXIS AND TWO-AXIS SOURCE 
SCATTERING EXPERIMENTS 

This chapter contains results of a series of scattering measurements. The target 
was a 61 cm X 6 cm X 2 cm brass plate. In the first series of measurements, the two-axis 
source was located 100 cm from the left side of the tank. The two-axis geophone was 
located 30 cm colinear and to the right of the source along a diameter of the tank (see Fig. 
4.1). The source was excited with a single cycle at 400 Hz. The phase difference between 
the two sources was 90°. 

First a background reading was recorded. Then the target was buried 85 cm from 
the source with its largest cross section facing the source-receiver. Its top edge was flush 
with the sand’s surface. So as viewed from the receiver the target was 61 cm wide and 6 
cm high. It took approximately 4 to 7 minutes to take the background signal, bury the 
target and then take the second reading. The results are presented in Fig. 4.2. The upper 
figure shows the vertical signal, the lower figure the horizontal signal. The vertical 
background signal shows two prominent features: one pulse lasting from approximately 
0.005 to 0.0235 and another from approximately 0.024 to 0.038 s. The signal with the 
target present was very similar. The only significant difference occurred at approximately 

^ 3 

0.02-0.027 s. Upon closer inspection, the two traces began to differ at about 0.015 s. 
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2-D Source Vertical and Horizontal 85 cm Signals 



Figure 4.2. Graph showing the output of the geophone with (dashed) and without (solid) 
the target present. The target range was 85 cm from the source foot. The two-axis 
source was driven with a single 400 Hz cycle 
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There was another, much weaker discrepancy at approximately 0.035-0.05s. This data 
therefore might indicate that the return from the target arrived at the receiver at about 
0.015s after transmission. The complete path distance was approximately 140 cm, so the 
propagation speed was approximately 90 m/s. The origin of the second pulse in the 
background and scattered signals is not known. 

The horizontal signal behaves very differently. The background and scattered 
signals are different right from the start. There is evidence of many more echoes or 
multipath arrivals than with the vertical signal. The horizontal signal does not lend itself to 
as simple analysis as the vertical signal. 

The measurements were repeated with the target 60 cm from the source. The 
signals are shown in Fig. 4.3. The vertical signal is similar to that in Fig. 4.2, except the 
discrepancy between background and scattered signals begins earlier, consistent with a 
shorter path length. The horizontal signal is also similar to that in Fig. 4.2. It is equally 
different than the vertical signal and equally difficult to interpret. 

The measurements were repeated with the one-axis source from Ref 1. The 
signals are shown in Figs. 4.4 and 4.5 for 85 cm and 60 cm target distances, respectively. 
The total time required to complete the two-axis and single-axis experiments was 
approximately 12 minutes. The background signals for the one and two-axis sources are 
roughly the same. The main difference was that the second pulse in the vertical signal was 
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Horizontal, 400 Hz, 60 cm 


2-D Source Vertical and Horizontal 60 cm Signals 




Figure 4.3. Same as Fig. 4.2 except for a 60 cm target range. 




Horizontal. 400 Hz, 85 cm 


1-D Source Vertical and Horizontal 85 cm Signals 
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Figure 4.4. Same as Fig. 4.2 except the source was the single axis source. 










1-D Source Vertical and Horizontal 60 cm Signals 










weaker with the one-axis source. It was surprising that the vertical signals with the target 
present were very similar. The identification of a feature that arrives earlier with the 60 
cm target than the 85 cm target is not obvious. Again the horizontal signals are different 
fi'om the vertical signals. Unfortunately, they are also difficult to interpret. 

The main conclusion fi'om these measurements is that is that the information in the 
vertical and horizontal signals was different. Conducting these measurements in a more 
fi*ee-field environment may yield more easily interpreted signals. Further, use of the 
information in both axis may perhaps lead to identification of surface wave components in 
the signals. 

Other methods used to extract the target fi’om the received obstructed signals 
included subtracting the obstructed and background signals. This process sometimes 
yielded a definite pulse that correlated to the buried object. However, most times, this 
process did not yield any useful information. This inconsistency was attributed to the 
theory that the propagation was dominated by reverberation. Ref 2 documents the effects 
of small changes in sediment particle size and the resulting interface that can make 
subsequent target detection difficult. It is not unreasonable to apply this theory to this test 
environment. Specifically, this test tank’s sides and bottom can be assumed to model rigid 
interfaces which can support many reflections and the resultant noise field can make 
isolation of the target reflection difficult. These reverberation effects become more 
evident when using analysis techniques designed to separate and identify reflections in 
acoustic propagation. This will be further described in the next section, Cepstrum 
Analysis. 
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B. CEPSTRUM ANALYSIS 


Autocorrelation and autospectrum analysis may identify the presence of two paths 
and the difference in their propagation times from the source to the geophone. When two 
or three paths are present, however, the identification of individual paths becomes more 
difficult [Ref 4], A Fourier transform of the output spectrum (called the cepstrum) will 
convert the individual interference components into more readily identifiable cepstral 
peaks at frequencies corresponding to Afyik = 1/1 Ti-Xk I, i ^ k = 1,2,3, 

..., r [Ref 4], Specifically, 


Cx = l/27cflog I X(e’“) I e*®" d® (4. 1) 

is the real cepstrum of a signal x. I X(e’’*) I is the magnitude spectrum of the signal x. 
Thus Cx is the inverse Fourier transform of the logarithm of the magnitude spectrum of 
signal X.. The limits of integration range from -n to n. 

The following example is taken from Ref 4 and demonstrates the 
potential of cepstrum analysis. A degraded sonar signal pulse 


t*exp-®‘ 


(4.2) 


47 



with damping coefficient a = 0.06, was degraded by adding a delayed, attenuated replica 
of the basic pulse. This data was input to the cepstrum algorithm. The results appear in 
Fig. 4.6. The upper trace in Fig. 4.6 shows the composite signal and the echo that was 
delayed 55 seconds relative to the basic waveform. This peak is easily seen in the lower 
trace of Fig. 4.6. Figures 4.7 through 4.10 show the results of cepstrum analysis applied 
to the scattering experiment data collected in section 4A. Each of these figures show the 
cepstrum of the background, target, and difference between the background and target 
signals. In general, the cepstrum plots show small amplitude features which do not 
correlate to any particular propagation path and the inherent difference between vertical 
and horizontal signals did not seem to have any effect on the resultant cepstrum signals. 

In each case, the cepstral signals followed the same form but did not show any clearly 
defined features that correlated to the target. The difference signals mirrored the target 
signals, but did not otherwise show any interesting features. The reasons for these results 
is not clear. They could be attributed to the reverberant field masking the potentially weak 
target feature (as described in Ref 2) or a modal feature of the tank that we do not 


understand. 
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Figure 4.7. Results of Cepstrum analysis of the vertical signal data shown in Fig. 4.3. The 
dashed line indicates the difference in the Cepstrum of the background and with the target 


Vertical Signal Results 
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Figure 4.8. Same as Fig. 4.7, except the data are from the horizontal geophone. 
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Horizontal Signal Results 






































Horizontal Signal Results 



Cepsrum For Target at 85 cm, 400 Hz Signal 
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V. CONCLUSIONS AND RECOMMENDATIONS 


This thesis describes the design and testing of a two-axis source and a two-axis 
receiver for use in buried object detection. The purpose of the two-axis source was to be 
able to selectively excite surface waves in water-saturated sand, exploiting the elliptical 
particle motion. The purpose of the two-axis receiver was to take advantage of the 
surface wave phasing to filter out signals fi-om non-surface wave modes of propagation. 

We successfully demonstrated the ability to generate particle motions with 
different phasing with the source. Definitive demonstration of surface wave generation 
was not possible, however, because of limitations imposed by due the size of the test tank 
Long pulses caused reverberation problems. Single cycle excitation resulted in different 
firequency and phase contents in vertical and horizontal signals, ruling out the possibility of 
detecting constant phase differences. 

Scattering experiments demonstrated the advantage of a two-axis receiver. The 
information content in the two channels were unexpectedly different. The main 
recommendation is to test the system on a beach to more closely approximate fi’ee-field 
conditions. 
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APPENDIX A. LABVIEW VI FOR DATA ACQUISITION 

One vi was used to acquire data. The vi is designed to be triggered by the HP 
3314A Function Generator. The vi should acquire a specific number of received 
waveforms fi'om each channel and average them. The three programming elements 
required were a “for loop,” shift registers and an intermediate vi call AI wave. The three 
elements can be seen in the simplified version shown in Fig. (A. 1). The book icon is a 
“for loop.” The arrows on each side are shift registers. In this case the “for loop” value is 
100. Therefore, 100 acquisitions are made, summed and then divided by the “for loop” 
value to get an average of the signal. This is timed and buffered acquisition. This means 
that the hardware clock is used to control the acquisition rate and the data is stored in an 
intermediate memory buffer after it is acquired from the input channel. Figure (A.2) 
shows the elements of the actual vi. The left side of Fig. (A.2) shows the front panel 
inputs controlling the number of iterations, number of channel inputs, specific device, 
trigger type, and scan rate information. Figure (A.3) shows the actual front panel. 



Figure (A. 1). Elements of Data Acquisition. 
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Figure (A. 3). LAB VIEW vi for data acquisition 
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APPENDIX C: FACT SHEET FOR SOURCE ELEMENTS 


WOOFER SPEAKER 
4"(10cm) Woofer 

FEATURES 

• Smooth frequency response 

• Light weight hard cone paper based on careful 
study ol internal loss values 

• Uprolied rubber edge minimises distortion and 
improves linearity 

• Long travel voice coil for high power handling 



SPECIFICATIONS 

Impedance 
Frequency response 
Free air resonance 
Ois 
vas 

Vas IS the volume of air 
equivalent to the compli 
SPL 

input power nominal) 
Input power (maximum) 
Magnet weight 
Speaker 


8 ohms 
50-7.000H2 
55H2+5H2 
0 35 
0 23 eft 

that IS acoustically 
iance of the cone 
84 + 2dB/W (im) 
5W 
tow 

228 g 8 0 02 
700 g 24 7 02 


HOW TO USE 

1 Installation: 



a Install the woofer unit from the enclosure front and keep it flush with the baffle surface 
b Secure the speaker so that no air leakage is permitted 

2. Connections-when using with a tweeter m a two-way system connect as shown 
Wnn\c< 



Figure (B. 1) Fact Sheet for geophones. 
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